The hunting for the discrete Painleve VI equation is over 
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Abstract 

We present the discrete, q-, form of the Painleve VI equation written as a three-point mapping and 
analyse the structure of its singularities. This discrete equation goes over to Pyi at the continuous 
limit and degenerates towards the discrete q-Py through coalescence. It possesses special solutions in 
terms of the g-hypergeometric function. It can bilinearised and, under the appropriate assumptions, 
ultradiscretised. A new discrete form for Py is also obtained which is of difference type, in contrast with 
the 'standard' form of the discrete Py- Finally, we present the 'asymmetric' form of q-Pyi as a system 
of two first-order mappings involving seven arbitrary parameters. 



The six Painleve equations were discovered a century ago by Painleve [1] who responded to the challenge of 
Picard [2] by classifying the integrable second-order equations of the form w" = f{w', w, t) and discovered 
the six transcendental equations that were subsequently named after him. These equations define new 
functions which extend into the nonlinear domain the special functions of the hypergeometric family, while 
constituting the nonautonomous extensions of the elliptic functions. Discrete analogues of these equations 
have been proposed much later, although, already at the time of Painleve, Laguerre [3] had examined 
integrable nonautonomous discrete systems and Shohat [4], half a century later, had proposed an equation 
which is known today as the discrete Pj. Discrete Painleve equations were subsequently derived in various 
contexts [5,6] and, finally, identified as such in [7]. All these equations were of difference type, where 
the independent variable n enters in an additive way. Multiplicative, g-discrete, Painleve equations were 
obtained for the first time in [8] where we presented the discrete, or g-discrete, forms of the first five 
Painleve equations. Subsequent work has led to the proposal [9] of a discrete analogue for all (but one) 
equations of the Painleve-Gambier [10] classification. 

Still, one element was missing: the discrete form of Painleve VI. A first step in this direction was 
accomplished when Jimbo and Sakai [11] discovered that the 'asymmetric', two-component, form of q- 
Pin has indeed Pyi as its continuous limit. Shortly afterwards the asymmetric forms of d-Piv and g-Py 
were also shown [12] to go over to Pyi at the continuous limit. (We must point out here that the latter 
discrete equation possesses more parameters than the asymmetric g-Pm and thus, as a discrete system, 
goes beyond Pyi)- Despite these interesting results, the question of the existence of a one-component, 
'symmetric' in the QRT terminology [13], form for the discrete Pyi remained open. In previous papers 
we have speculated on the form of this equation [14] and its singularity structure [9], but its precise form 
was missing. This is remedied in this paper. 

We shall start by presenting the form of the discrete Pyi which is of multiplicative, q-, type. Consequently 
we shall refer to it as g-Pyi. Once the form is given we shall verify that it has all the required properties: 
a continuous limit to Pyi, special solutions of q- hypergeometric type, a degeneration through coalescence 
to g-Py and, foremost, a generic, symmetric, singularity pattern (as expected from our work in [9]). The 
bilinearisation of g-Pyi, its generalisation to as asymmetric form and its ultradiscretisation [15] complete 
this paper. One caveat is in order at this point. No systematic derivation of q-Pyi will be presented. Its 
form was not derived from first principles but, rather, through intuition and inspiration, two invaluable 
investigation ingredients. 

Let us give the form of q-Pyi- With the usual notation x = x(n), x = x(n + 1), x = x(n — 1), we have: 

(xx — zz) (xx — zz) {x — az){x — z / 'a){x — bz)(x — z /b) 
(xx-l)(xx-l) ~ (x -c)(x- l/c)(x -d)(x- ijd) ^ ' 

where z = z X n (and, of course, z = Xz, z= z/X) and a,b,c,d are free constants. The singularity 
pattern of equation (1) can be obtained in a straightforward way. First, suppose that at a given iteration 
we have x — c while x has a generic value. This means that the rhs of (1) diverges. This has as a 
consequence xx = 1 and thus x = 1/c (i.e. a vanishing denominator at the next step), leading to regular 
values for x thereafter. Similarly, if x goes through a zero of the numerator, x — az for instance, we 
must have xx — zz and thus x = z/a. Thus the numerator will vanish at the next step ensuring the 
confinement [16] of this singularity. Finally, from the form of (1), it is clear that x can go through the 
value zero. A detailed analysis of this situation leads to the condition zz_ — z 2 , in which case x — is 



not a singularity at all. Thus we find z = z A". The net result of this analysis is that the singularities 
of g-Pyi are grouped in pairs leading to 8 singularity patterns: {az,z/a} 7 {z/a, az} 7 {bz, z/b}, {z/b, bz}, 
{c, 1/c}, {1/c, c}, {d, 1/d}, {1/d, d} where confinement occurs in just one step. This is precisely what we 
have postulated in [9] where we stated that q-Pyi, just as its continuous counterpart should possess a 
symmetric and generic singularity pattern. 

Equation (1) is indeed a discrete form of Painleve VI as can be assessed from the continuous limit. 
Putting A = e e , a — —e ea , b — e e/3 , c = — e C7 , d — e eS we find in the limit e — > 0: 

d 2 x 1 / 1 1 1 1 \ ( dx\ 2 (1 1 1 1 I \ dx 

+ 



dz 2 2 \x + 1 x — 1 x + z x — z J \dz J \z z — 1 z + 1 x — z x + z ) dz 

(x 2 z 2 )(x 2 - 1) ( {a 2 -l/A)z Q3 2 -l/4)z 7 2 m 
z 2 (z 2 -l) v (x + z) 2 (.x-z) 2 (x + l) 2 (x - l) 2 ) {> 

which is indeed Painleve VI in non canonical form. In order to make it canonical we introduce the 

following change of variables: z = (1 + v / C)/( 1 ~ VC) an d x — (VC + W )/(V( — w ) an d obtain: 

d 2 w 1/1 1 1 \ f dw\ 2 (\ 1 1 \dui 



+ 7 + > - t + ^t + 



w(w _i )(w _ c) / 2 _ , ^ 2 (C-i) , (i-/3 2 )C(C-i) A 
2C 2 (C-i) 2 V w 2 + ( w -i) 2+ (^-C) 2 /' 

The discrete Pyi equation falls nicely into the pattern of degeneration through coalescence that has been 
established for the continuous Painleve equations and which has been verified for the discrete Painleve as 
well [8,17]. Indeed, the proper limit of g-Pvi allows us to recover q-Py. Following our customary notation 
[18], we shall denote by uppercase letters the variables of g-Pyi and by lowercase those of q-Py. We put 
X = x, Z = z/S, A = a/5, B = b/5, C = c, D = d and then take the limit 5^0. We obtain thus: 

(x - c)(x - l/c)(x - d)(x - 1/d) . . 

(xx-1 ){xx-l =" 77 7171 t —, s 4 

(1 - ax/z)(l — bx/z) 

which is precisely q-Py in canonical form. 

The continuous Pyi equation has solutions in terms of hypergeometric functions for some special values 
of the parameters. The same holds true for g-Pyi- The simplest way to obtain these special solutions 
is to use the splitting technique [19]. We separate equation (1) in two discrete Riccati, homographic, 
equations in the following way: 

xx — z~z (x — az)(x — bz) . . 

(5a) 



xx — 1 (x — c)(x — d) 
xx — zz (x — z/a)(x — z/b) 



(5b) 



xx — 1 (x — 1/c) (a; — 1/d) 
The two equations of system (5) are indeed homographic and compatible provided the condition ab = Xcd 
holds. The linearisation of the discrete Riccati is obtained through a Cole-Hopf transformation x = P/Q, 
resulting to the linear equation 

Q(az — d)(az — c)((a + b)z — c— d) 

+aQ ((a + b)z((cd\ - l)z 2 + \-cd)-(c + d)((ab - l/X)z 2 + 1 - ab/X)) 

-Q(a-dz)(a-cz)((a + b)z-c-d) = (6) 

Equation (6) has the hypergeometric equation as continuous limit. This limit is simpler to obtain if we 
start from the discrete Riccati (5) and implement the continuous limit introduced above. With the same 



ansatze we are thus led to a continuous Riccati equation which can be linearized through the Cole-Hopf 



u:i -C)0 + (/3-<5-(/3 + 7 + i)C)^- hu =0 (7) 



i.e. precisely the Gauss-hypergeometric equation. Just like continuous Pyi, equation (1) has also rational 
solutions. The simplest one, x — y/z, is obtained provided the conditions a — c and b — d hold. Higher 
solutions of both linearisable and rational type can be obtained throught the application of the Schlesinger 
transformations of g-Pyi, which are currently under study [20]. 

Another interesting aspect of g-Pvi is its bilinearisation. In [14] we have explained how the bilinear, 
Hirota, form can be obtained using the singularity pattern as a guide. Indeed, we expect the number of 
r-functions to be equal to the different singularity patterns [21]. We start by introducing 8 r-functions 
t oj n, t c , T d , a ai (J;,, cr c , a d and regroup them in combinations <pi = T a a a , <f> 2 — f a a a , 4> 3 = T b a b , 
etc. Parallely we introduce the 8 quantities oti = z/a, a 2 = za, a 3 = z/b, etc. The bilinearisation is 
introduced by the following ansatz for x: x = cti — 4>i/ip, i = 1, . . . , 8 where ip is, a priori, bilinear in the 
r-functions r, a. From the 8 possible definitions of x we obtain, after elimination of ip, 6 independent 
bilinear equations. We can write them in a most symmetric form as: 

<t>i{aj - a k ) + (pj(a k - a % ) + 4>k{ai - «j) = (8) 

where we have used the fact that tjj = (<fii — 4>j)/(ai — aj) for all i ^ j. There are 56 such equations which 
can indeed be obtained from 6 independent ones, for instance those with i = 1, j = 2 and k = 3, . . . , 8. 
The remaining two bilinear equations can be obtained by substituting the ansatz for x into (1). The net 
result, after all possible simplifications, is an equation that can be separated into two bilinear ones: 

aa a T a -^g_ a T a = K{a-^-)<j b T b (9a) 

aCT c r c - ^ct c t c = K (c - ^)a d r d (96) 

Equation (9a) can be written in a different way where a and b are permuted, but this two forms arc 
equivalent, provided equations (8) hold. The same holds true for equation (9b) under a permutation of 
c and d. It is clear from the form of equations (9) that the value of K can be modified through a gauge 
transformation which leaves (8) invariant. Thus, with the appropriate choice of gauge, we can take K = 1 
without loss of generality. The system that consists of equations (9) and 6 independent equations among 
those of (8) is the bilinear form of g-Pvi- 

One final point concerning g-Pyi is its ultradiscretisation. This method has been introduced by the 
Tokyo-Kyoto group [15] as a method for obtaining discrete equations where, provided one ensures the 
right parameters and initial conditions, the dependent variable takes only integer values. Thus the 
evolution introduces a generalised cellular automaton. The ultradiscretisation of g-Pyi is obtained using 
the standard techniques of [22]. We first rewrite (1) as: 

_ (x + az)(x + z/a)(x + bz)(x + z/b) 
V - ~ (x + c)(x + l/c)(x + d)(x + 1/d) ( ' 



by inverting the sign of the constants a, 6, c, d. Next, we introduce A = e 1 / 6 , x = e x / e , y = e y / £ and similar 
ansatzc for z, a, b, c, d. The limit of (10) which uses extensively the identity lim e „» log(e a / c + e' 3 / 6 ) = 
max(a,/3), leads to the ud-Pyi: 

Y + Y = ma,x(X, Z + A) + max(X, Z - A) + max(X, Z + B) + max(X, Z - B) 

- max(X, C) - max(X, -C) - max(X, D) - max(X, —D) (11a) 

X + X = max(r,2Z+ 1) - max(Y", 0) (116) 

The study of Pyi leads naturally to more equations that we will briefly examine here. The first is a 
'symmetric', one component, difference equation which can be obtained almost by inspection from g-Pyi, 
by considering a different coalescence. Indeed taking X = 1 + 5x, A = 1 + 5k, A = 1 + 5a etc., we find: 

(x + x — z — z)(x + x — z — z) (x — z — a)(x — z + a)(x — z — b)(x — z + b) 
(x + x)(x + x) (x — c)(x + c)(x — d)(x + d) 

where z = z a + nn. This equation is a new discrete form of Py- This can be assessed from the continuous 
limit of (12). Putting z = VC, x = \fC/(l - w) and a = ae 2 , b = j3/e 2 , c = -ye 2 , d = f3/e 2 + 5 we obtain: 

d 2 w ( 1 1 \ / dw\ 2 1 dw (w — l) 2 7 5w w(w — 1) 

^C 2 = \2w~ + w~l) \d() _ C^C + 2C 2 w' + 2f3 2 {w-iy (13) 

Finally we can ask whether there exist asymmetric (two-component) extensions of both (1) and (12). It 
turns out that this is straightforward. We obtain thus the asymmetric g-Pyi 

(yx- zz)(xy - zz) (y - az)(y - bz){y - cz){y - dz) 



(yx- l)(xy- 1) {y - p)(y - q)(y - r)(y - s) 

(xy — zz)(xy — zz) (x — z/a)(x — z/b)(x — z/c)(x — zjd) 



(14a) 
(146) 



(xy - l)(xy - 1) (x - l/p) (x - l/q) (x - 1/r) (x - 1/s) 

where z = zV\, z = z/\f\, and with the constraints abed — 1 and pqrs = 1. Similarly, we find the 
asymmetric d-Py 

(y + x — z — z) (x + y — z — z) (y — z — a)(y — z — b)(y — z — c)(y — z — d) 



(y + x)(x + y) (y - p)(y - q)(y - r)(y - s) 

(x + y - z - z)(x + y - z - z) (x - z, + a)(x - z, + b)(x - z. + c)(x - z + d) 



(15a) 



(156) 



(x + y)(x + y) (x +p)(x + q)(x + r)(x + s) 

where z = z + re/2, z = z — k/2, and with the constraints a + b + c + d = and p + q + r + s = 0. These 
equations are the most general discrete Painleve equations known to date. They have seven discrete 
parameters and thus they lie beyond the already known asymmetric q-Pm, d-Piy and g-Py- Still, we 
expect their continuous limits to go over to Pyi- The fact that these equations do not have the constraints 
of the symmetrical form make us believe that they possess the property of self-duality just like asymmetric 
g-P ln [23] and asymmetric d-Ppv and q-P v [24]. Thus it is interesting to study these equations in the 
framework of what we call the Grand Scheme i.e. the simultaneous analysis of a discrete equation and 
its Schlesinger transformations. This aspect of q-Pvi and d-Py is the object of a future study [20]. 
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